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We present the first direct measurement of the top-quark mass using tt events decaying in the 
hadronic r + jets decay channel. Using data corresponding to an integrated luminosity of 2.2 fb _1 
collected by the CDF II detector in pp collisions at ^/s — 1.96 TeV at the Fermilab Tevatron, 
we measure the tt cross section, <r t j, and the top-quark mass, M top . We extract M top from a 
likelihood based on per-event probabilities calculated with leading-order signal and background 
matrix elements. We measure a t f = 8.8 ± 3.3 (stat) ± 2.2 (syst) pb and Mt op = 172.7 ± 9.3 (stat) ± 
3.7 (syst) GeV/c 2 . 

PACS numbers: 14.65.Ha, 13.35.Dx 



The mass of the top quark, M top , and the top-quark 
pair production cross section, <7 t j, have been extensively 
studied at both the Fermilab Tevatron and the Large 
Hadron Collider at CERN [1 3 . However, final states 
of the top-quark decay that include a tau lepton (r) are 
relatively unexplored due to the difficulty of identifying 
the tau and rejecting quantum chromodynamic (QCD) 
processes that can mimic its hadronic decay mode. The 
top quark and the tau belong to the heaviest third gen- 
eration of standard model (SM) fermions and may play a 
special role in electroweak symmetry breaking. Discrep- 
ancies from the SM expectation in these unexplored de- 
cay channels could point to new top-quark sector physics. 

In this Letter, we present the first direct measurement 
of the top-quark mass in the hadronic r + jets decay 
channel (r + jets) |4]. We measure a t t an d M top us- 
ing data corresponding to 2.2 fb _1 of integrated lumi- 
nosity collected by the CDF II detector [5] in pp colli- 
sions at = 1.96 TeV at the Fermilab Tevatron. The 
DO Collaboration previously measured a t i in the r + 
jets decay channel with data corresponding to 1 fb _1 of 
integrated luminosity to be 6.9 ± 1.5 pb [BJ, assuming 
-Mtop = 170 GeV/c 2 , and reached a signal purity of 52%. 

Assuming three generations in the SM and a unitary 



quark-mixing matrix, the top quark decays almost ex- 
clusively to a W boson and b quark. We select pair- 
produced top-quark events in which one of the W bosons 
decays into a pair of light quarks and the other decays 
to a tau and a neutrino. This decay channel represents 
15.2% of the tt branching ratio and results in a final state 
with a tau, a neutrino, two b quarks, and two light-flavor 
quarks (u, d, or s). Although the tau can decay lep- 
tonically to an electron (e) or muon (/it) and a pair of 
neutrinos, these events are difficult to differentiate from 
electrons or muons from W boson decays. As a result, 
we select events with the tau decaying to a neutrino and 
a narrow jet of hadrons, which are usually charged and 
neutral pions, that correspond to 9.8% of all tt decays. 
We use an artificial neural network (NN) to reduce the 
QCD multijet background contribution. The additional 
neutrino produced in the tau decay complicates the tau 
reconstruction. To solve this, we adapt a missing mass 
calculator method [7] to the r + jets topology to infer 
a unique solution for the neutrino four-momentum with 
sufficient precision to reasonably reconstruct the tau. We 
use a binned likelihood fit based on the predicted and ob- 
served number of events to measure a t i . Then, to extract 
■Mtop j we use a likelihood function built from signal and 
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background probabilities calculated with the predicted 
differential cross sections for tt and W + four-parton pro- 
duction, respectively. 

The data used in this measurement are selected using 
a multijet online selection (trigger) that requires at least 
four calorimeter energy clusters with transverse energy 
[5] (Et) greater than 15 GeV each and a total sum Et 
of all clusters greater than 175 GeV [§]. Jets are recon- 
structed by a cone algorithm that clusters energies in 
calorimeter towers within a fixed cone size of AR = 0.4 
[TP] where AR = yj Arj 2 + A(f> 2 0. In the offline anal- 
ysis, events are required to have exactly four jets with 
Et > 20 GeV, missing transverse energy (^t) greater 
than 20 GeV, and a single hadronically-decaying tau se- 
lected as described below. Jet energies are corrected for 
nonlinearity of the detector response and multiple pp in- 
teractions within the bunch crossing [TT]. One of the 
four jets must be identified as having originated from a 
b quark (6-tagged) using a secondary vertex finding algo- 
rithm 12 . Hadronically-decaying taus appear as narrow 
jets with an odd number of tracks and low neutral pion 
multiplicity. We select taus using a two-cones algorithm 
[T5] . The inner cone defining the signal region has a size 
set to the lesser of 10° (0.17 rad) and (5 GeV)/.E c i rad, 
where E c \ is the energy of the calorimeter energy cluster 
associated with the tau candidate. The second cone with 
a size of 30° defines an isolation region outside of the 
signal cone. A tau must have one or three tracks in the 
signal region and no tracks in the isolation region. We 
require the Et of the tau energy cluster to exceed 20 GeV 
and the Et of the visible tau to exceed 25 GeV, where 
visible refers to the combination of the track and neutral 
pion information. We require that the calorimeter energy 
in the isolation region be less than 10% of the visible tau 
energy. Finally, we veto events with an identified electron 
or muon. 

The dominant background for this analysis is high jet 
multiplicity QCD events with a jet misidentified as a tau. 
To reduce this background, we develop a NN to distin- 
guish between tt — > r + jets and QCD multijet events. 
The NN is trained using QCD multijet events, obtained 
from data by selecting events with a tau candidate with 
at least one track in the isolation region and passing all 
other selection requirements, and tt events generated us- 
ing the pythia Monte Carlo generator [TT] coupled with 
a GEANT [TH] based CDF II detector simulation [IB] . To 
properly account for tau polarization effects, the tau de- 
cays are modeled by the tauola package [17J . The NN 
uses eight variables that exploit the topological differ- 
ences between QCD multijet and tt events including Et 
and the sum of the Et of various combinations of the tau 
and jets [4]. After training the NN, we find good separa- 
tion between QCD multijet and tt events. Optimal sig- 
nal significance, defined as the number of expected signal 
events divided by the square root of the total number of 
observed events, is achieved by removing events associ- 



ated to a NN value below 0.85. We initially select 162 
events of which 41 events survive the 0.85 NN require- 
ment. 

Due to the difficulty in simulating QCD multijet 
events, b quark tagging algorithms, and the production of 
heavy flavor quarks in association with W bosons, we es- 
timate the background contributions with a data-driven 
approach similar to that described in Ref. [JJ]. We use 
the NN output distribution to fit the contributions of 
the signal and background processes to the data. This 
is done both before and after applying 6-tagging require- 
ments. Since most of the selected data events return 
a NN value below 0.7, the fit is dominated by events 
outside of the signal region. We begin by calculating 
the contributions of the signal and each background pro- 
cess before the 6-tagging requirement is applied. The 
tt and electroweak background contributions are deter- 
mined from simulation. Diboson, single top-quark, and Z 
+ jets production are modeled using pythia, madevent 
[18] , and ALPGEN [19] respectively, with PYTHIA used for 
parton showering and underlying event generation. Each 
of these processes' contributions is set to its expectation 
based on its respective theoretical cross section [30] [2JJ , 
the total integrated luminosity of the data, and the accep- 
tance determined from simulation. The tt contribution 
is modeled with pythia and is similarly normalized us- 
ing the next-to-next-to-leading order SM tt cross section 
prediction [52]. For all simulated events, a GEANT based 
simulation is used to model the CDF II detector response. 
With these contributions known, we determine the con- 
tributions from QCD multijet and W + jets events by 
fitting the shape of the NN output distribution for each 
component (with the previously calculated contributions 
fixed) to the data before applying the NN selection and 
6-tagging requirements. The QCD multijet sample is se- 
lected from data as previously described while the W 
+ jets events are modeled with ALPGEN similarly to Z 
+ jets events. We fit these distributions with a binned 
Poisson likelihood as seen in Fig. [I] The contribution 
of QCD multijet events in the signal region is calculated 
from this fit. All remaining events arc assumed to come 
from W + jets production. 

For each process except QCD multijet events, the con- 
tribution after applying the ^-tagging requirement is cal- 
culated by applying 6-tagging efficiencies measured in 
Ref. [T5] to the initially calculated contribution. Incor- 
rect tagging of light quarks and tagging of the b and c 
quarks have inherently different uncertainties. To prop- 
erly estimate uncertainties associated with the contribu- 
tion of W + jets events, this contribution is divided into 
W + light flavor (W+lf) and W + heavy flavor (Wbb, 
Wcc, and Wc) parts with separately-estimated uncer- 
tainties. To calculate the contribution from QCD mul- 
tijet events, we apply the ^-tagging requirement to the 
QCD multijet sample. We then combine tt and the other 
background processes into a single sample with the rel- 
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TABLE I. Predicted number of selected r + jets candidate 
events from each considered process after all selection re- 



quirements are applied assuming cr t j = 7.4 pb and Af top = 
172.5 GeV/c 2 . The uncertainty is a combination of statistical 
and systematic uncertainty. 

Source Number of events 

Diboson 0.19 ±0.01 

Single top quark 0. 16 ± 0.01 

Zbb 0.29 ± 0.04 

Wbb 0.6 ± 0.5 

Wcc 0.3 ±0.3 

Wc 0.2 ±0.1 

W±lf 0.5 ±0.6 

QCD multijet 18.2 ±4.1 

Total bkgd 20.4 ±4.2 

it 18.2 ±2.8 

Total predicted 38.6 ± 5.0 

Observed 41 



ative contributions fixed to their calculated values. The 
NN output distributions of these two samples are then 
fit to the data selected with the 6-tagging requirement, 
and the contribution of QCD multijet events in the signal 
region is derived from the result. Each process's contri- 
bution, assuming a t i = 7.4 pb and M top = 172.5 GeV/c 2 , 
is given in Table [TJ Of the 41 selected events, we expect 
roughly 18 QCD multijet events and 18 ti events. From 
simulation studies, we estimate that 76.5 ± 0.5% of the 
selected ti events correspond to a hadronic r ± jets final 
state. The other major contributions to the ti events are 
all-hadronic ti decays (12.3 ± 0.4%) and ti — > e + jets 
(5.3 ± 0.3%). 




-0.2 0.2 0.4 0.6 0.8 1 1.2 
NN output 



FIG. 1. Fit to the NN output shape before applying the b- 
tagging requirement. The arrow marks the lower bound on 
the signal region. 

We measure a t i using a likelihood function based on 
a Poisson probability distribution comparing the num- 
ber of observed (N a ) and predicted (N p ) events for a 



given a t i written as L = e~ Np Np°/(N \). We consider 
the negative logarithm of this function over values of 
a t i from 5 to 15 pb where N p is recalculated at each 
point with the fraction of QCD multijet events kept con- 
stant to the value calculated for er t j = 7.4 pb. The re- 
sult is fit with a 2nd order polynomial which is mini- 
mized to extract the central value and statistical uncer- 
tainty. The cross section value determined by the fit is 
a tI = 8.8 ± 3.3 (stat) pb. 

The dominant sources of systematic uncertainty in- 
clude the acceptance, selection efficiencies, background 
estimate, and luminosity. For acceptance effects, we con- 
sider uncertainties on the jet energy scale (JES) [TT] (0.6 
pb), parton showering models (0.5 pb), parton distri- 
bution functions (PDF) (0.5 pb), initial and final state 
radiation (ISR, FSR) (0.5 pb), and color reconnection 
[2"31 (0.4 pb). We consider systematic uncertainties on 
the efficiency measurements from the 6-tagging (0.4 pb), 
tau identification (0.2 pb), and trigger efficiency (0.1 
pb) scale factors. The background systematics come 
from the W + heavy flavor scale factor uncertainty [12] 
(0.1 pb) and the QCD multijet contribution, which is 
the dominant systematic uncertainty. We measure this 
uncertainty (1.8 pb) by comparing the NN output dis- 
tribution shapes of the QCD multijet events and data 
dominated by QCD multijet events which are selected 
by removing the ^ requirement. Finally, the uncer- 
tainty on the integrated luminosity is 6% [M] (0.5 pb). 
Combining all these sources in quadrature results in the 
total systematic uncertainty of 2.2 pb, a 25% uncer- 
tainty. We measure a t i assuming M top = 172.5 GeV/c 2 
to be 8.8 ± 3.3 (stat) ± 2.2 (syst) pb, which is consistent 
with the next-to-next-to-leading order SM prediction of 

7A5±° f 3 P b m- 

We calculate M top from a likelihood function based on 
probabilities corresponding to the signal and background 
hypothesis for each event. These probabilities are calcu- 
lated from the differential cross section for ti and W + 
four-parton production, respectively. The method uses a 
similar approach to the previous measurement in the e 
and li + jets decay channels [25]. The signal probability 
is based on a ti leading-order matrix element which as- 
sumes qq production |26j and is calculated over 31 input 
mass values ranging from 145 to 205 GeV/c . Since it 
does not depend on M top , the background probability is 
calculated once for each event using a W + four-parton 
matrix element from the VECBOS [27] generator. 

The tau decay adds an extra complication by introduc- 
ing a second neutrino in the event. We reconstruct this 
additional neutrino by adapting a method developed for 
the reconstruction of a resonance decaying to tt [7J to 
the r ± jets topology. We find from simulation studies 
of ti — » r ± jets events that the neutrino from the tau 
decay is nearly collinear with the hadronic components 
of the tau decay as their 9 and <j> angles tend to agree 
within 0.1 radians. Additionally, we find that the <fi an- 
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gle of the neutrino from the W boson is within 1 radian 
of the 4> direction associated with the $t- Simulation 
studies show that these statements hold true in greater 
than 99% of events. We introduce a four-dimensional 
scan over the angles of both neutrinos restricting them 
to the above ranges. Assuming the neutrino mass is neg- 
ligible and that the W boson and tau have masses of 
80.4 GeV/c 2 and 1.777 GeV/c 2 , respectively, we com- 
pletely solve for the four- vector of each neutrino for each 
set of angles in the scan. We then compare the pre- 
dicted x- and y-components of the |5t from the neutrino 
solutions to the measured $t components with Gaus- 
sian probability functions and choose the set of angles 
that returns the greatest probability. This method accu- 
rately reconstructs the four-momentum of the neutrino 
from the tau decay, but it does not perform as well with 
the four-momentum of the neutrino from the W boson 
decay. Therefore, we use the result of this method only 
to determine the tau four-momentum in the event, while 
the neutrino from the W boson decay is reconstructed in 
the method as it would be for the e or /i + jets channel 
[25j by assuming the tt system is produced with no pt- 
Studies with simulated electron + jets events in which 
the electron energy is smeared to match the tau energy 
resolution show that any bias introduced by these meth- 
ods is removed by the final calibration 0]. 

Each probability is calculated by integrating over the 
differential cross section for the appropriate process: 

P=IJ da(y)f(Qi)f(Q~2)W(x, y)dq x dq 2 , (1) 

where da is the differential cross section, / is the parton 
distribution function (PDF) for a quark with momentum 
fraction of the incident proton q, x refers to observed 
quantities, y refers to parton level quantities, and W(x, y) 
is the transfer function used to map x to y based on sim- 
ulation studies. The event probability is a sum of the 
signal and background probabilities weighted by the sig- 
nal and background fractions, respectively. To improve 
the statistical uncertainty on the M top measurement, the 
likelihood function includes a Gaussian constraint on the 
background fraction set to 0.5 ±0.1 from Table [TJ We 
evaluate the likelihood function for each of the 31 input 
top-quark masses and fit the result with a second-order 
polynomial to derive Af top and its statistical uncertainty. 
We calibrate the measurement on 21 simulated tt samples 
covering a mass range of 155 to 195 GeV/c 2 . The like- 
lihood function and fit for the data before applying the 
calibration functions can be seen in Fig. [2] We measure 
M top to be 172.7 ±9.3 (stat) GeV/c 2 . 

The largest systematic uncertainty comes from the 
JES and is calculated to be 3.4 GeV/c 2 . We also con- 
sider systematic uncertainties from the differences in par- 
ton showering models (0.5 GeV/c 2 ), color reconnection 
(0.5 GeV/c 2 ), ISR and FSR (0.3 GeV/c 2 ), PDF's (0.1 
GeV/c 2 ), and the uncertainty on the fraction of tt pairs 




150 160 170 180 190 200 
M lop (GeV/c 2 ) 

FIG. 2. Negative log of the top-quark mass likelihood as a 
function of Af top for all data events. The calibration functions 
have not yet been applied. 



produced from gg fusion (0.2 GeV/c ). The background 
fraction uncertainty is measured by shifting each back- 
ground source within its uncertainty from Table [I] (0.5 
GeV/c 2 ). We consider uncertainties from different 6-jet 
fragmentation models and semileptonic branching ratios 
for jets from b quarks as well as shifts in the energy scale 
of these jets (0.4 GeV/c 2 ). We also account for shifts 
from the tau energy scale (0.2 GeV/c ). The pileup sys- 
tematic uncertainty (1.0 GeV/c 2 ) accounts for a known 
mismodeling in the luminosity profile of the simulation. 
Uncertainty due to local non-linearity of the method and 
any assumptions used is estimated by shifting the cali- 
bration function within its uncertainty (0.2 GeV/c 2 ). We 
take the remaining 0.14 GeV/c uncertainty on the fit of 
the mass residual (defined as the true mass subtracted 
from the measured mass) across all 21 mass points as an 
uncertainty on the limited size of the simulation sample. 
We find Aftop to be 172.7±9.3 (stat) ±3.7 (syst) GeV/c 2 
in agreement with the most recent Tevatron combination 
of 173.18 ±0.94 GeV/c 2 [3]. 

Using data corresponding to an integrated luminosity 
of 2.2 fb _1 , we have made the first direct measurement 
of the top-quark mass in tt events identified as decaying 
to a hadronic r + jets topology. Assuming a top-quark 
mass of 172.5 GeV/c 2 , we find the tt pair production 
cross section to be 8.8 ± 3.3 (stat) ± 2.2 (syst + lumi) 
pb. This value is consistent with the next-to-next-to- 
lcading order SM prediction [35] and recent measure- 
ments [28 , including the DO measurement in the same 
decay channel [5] . We measure the top-quark mass to be 
172.7 ± 9.3 (stat) ± 3.7 (syst) GeV/c 2 in agreement with 
the Summer 2011 Tevatron top-quark mass combination 
of 173.18±0.94 GeV/c 2 [3j. These measurements demon- 
strate that we can do complex analyses with tau leptons 
even in a high jet multiplicity environment at hadron col- 
liders. This is particularly interesting at the LHC where 
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new physics, e.g. SUSY, could preferentially lead to final 
states with tau leptons. 
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